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ABSTRACT : The experiment involved observations of 2,514 Holstein-Friesian cows to determine the effects of environmental 
factors (cow's age, calving season, weight and sex of calves, housing system) and genetic factors on gestation length in dairy cattle and 
the correlation between gestation length and other reproductive traits (calving ease, stillbirth rates and placental expulsion). Genetic 
parameters were estimated based on the sires of calved cows (indirect effect) and the sires of live-born calves (direct effect). The 
following factors were found to contribute to prolonged gestation: increasing cow's age, male fetuses and growing fetus weight. Optimal 
gestation length was determined in the range of 275-277 days based on calving ease and stillbirth rates. The heritability of gestation 
length was estimated at 0.201-0.210 by the direct effect and 0.055-0.073 by the indirect effect. The resulting genetic correlations suggest 
that the efforts to optimize (prolong) gestation length could exert an adverse influence on the breeding value of bulls by increasing 
perinatal mortality and calving difficulty. The standard errors of the investigated parameters were relatively high, suggesting that any 
attempts to modify gestation length for the purpose of improving calving ease and reducing stillbirth rates should be introduced with 
great caution. (Key Words : Holstein-Friesian Breed, Calving Ease, Stillbirth Rates, Heritability, Genetic Correlations) 



INTRODUCTION 

Gestation length (GL), the period from effective 
fertilization until calving, is a reproductive trait that 
significantly affects cattle breeding and production. The 
cow's age is the key environmental factor influencing 
gestation length. Gestation length is shorter in heifers than 
in older cows (Przysucha and Grodzki, 2009). High 
temperatures in the summer speed up calving and shorten 
gestation length (McClintock et al., 2003). High milk yield 
prolongs gestation, as suggested by a positive genetic 
correlation between GL and milk production levels (Silva et 
al. 1992). Both longer and shorter gestation periods 
contribute to a higher number of stillbirths (Norman et al., 
2009). GL values should be analyzed in conjunction with 
calving ease (CE), stillbirth rates (SB) and placental 
expulsion (PE) (Hansen et al., 2004; Przysucha and Grodzki, 
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2009). Estimated heritability of GL based on the direct 
effect were determined in the range of h 2 = 0.27-0.45, 
whereas much lower values were reported based on the 
indirect effect (Hansen et al., 2004; Jamrozik et al., 2005; 
Norman et al. 2009). The above suggests that the main 
genetic factors affecting GL heritability are the fetus and, to 
a smaller extent, it's mother (Nogalski, 2003). The fetus' 
dominant effect on GL could be related to the induction of 
parturition. The fetal adrenal cortex secretes Cortisol, 
increasing placental prostaglandin levels which initiate 
parturition (Stabenfeldt et al., 1992). Nadarajah et al. (1989) 
and Cervantes et al. (2009) reported significant genetic 
correlations between GL vs. CE and SB, whereas Hansen et 
al. (2004) observed a weak correlation between GL and 
other reproductive traits. Owing to the high significance of 
CE and SB values in the production of high-yielding dairy 
cattle, other traits contributing to the optimization of 
reproductive traits should also be studied. 

The objective of this study was to determine the effect 
of selected factors on gestation length in dairy cattle and the 
correlations between gestation length and other 
reproductive traits. 
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MATERIALS AND METHODS 

The study investigated 2,514 Black-and-White Polish 
Holstein-Friesian cows kept in four cattle sheds, including 
two sheds in a tied-stall system and two sheds in a free-stall 
system. Gestation length (GL) was determined as the 
number of days between effective fertilization until calving. 
The animals were divided into five classes based on the 
reported GL values: up to 275 days, 275 to 277 days, 278 to 
280 days, 281 to 283 days, and above 283 days. Based on 
the age factor, the cows were divided into primiparous, 
second-calf and third-calf females. Three calving seasons 
were identified: fall-winter (1 October to 31 January), 
winter-spring (1 February to 31 May) and summer (1 June 
to 30 September). The variability in calving ease was 
determined based on direct observations of calving ease 
(CE), stillbirth rates (SB) and placental expulsion (PE). 
Twin pregnancies were not taken into account. Calving ease 
was described in four categories: 1 - unassisted calving, 2 - 
calving with little intervention (assisted by one person), 3 - 
calving with considerable intervention (assisted by several 
persons), 4 - very difficult calving (hard pull) requiring 
veterinary intervention. Stillbirth values were assigned the 
following scores: 1 - born alive, survived the first 24 h after 
birth, 2 - dead at birth or died within the first 24 h after birth. 
Placental expulsion was classified as follows: 1 - 
spontaneous expulsion within 12 h after birth, 2 - 
spontaneous expulsion within 24 h after birth, 3 - placenta 
removed by a veterinarian. 

Distribution (%) of calving ease, stillbirths and 
placental expulsion in respect of investigated factors is 
presented in Table 2. Calf sex was described as: 1 - female 
calf, 2 - male calf. Calves were weighed directly after birth 
with the accuracy of 0.1 kg. They were divided into four 
weight groups: up to 34 kg, 35 to 37 kg, 38 to 40 kg, and 
above 40 kg. 

The effect of GL on CE was evaluated using the x test. 
We wanted to examine whether the length of pregnancy has 
an impact on the share of dystocia in dairy cattle. 



X 



fi-Fi)2 
Fi 



2008). A general linear model was developed in view of the 
key investigated traits as well as first-order interactions. 
The following classification model was used: 

yijklmn 

= n+Hi+Wj+S k +Ui+P m 
+(HxW)ij+ .. +(UxP) lm +e Ijk]mn 

where: 

Yijkimn = phenotypic value of trait, 

|i = general average, 

H, = fixed herd effect (i = 1, .., 4), 

Wj = effect of successive calving (j = 1, .., 3), 

St = effect of calving season (k = 1, .., 3), 

Uj = effect of housing system (1=1, 2), 

P m = effect of calf sex (m = 1, 2), 

(HxW)^, .., (UxP) lm , = interactions between factors, 



~ijklmn ' 



: random error. 



where: 

fi = numbers obtained 
Fi = expected number 

The observed CE, SB and PE values were not 
characterized by normal distribution, therefore, they were 
subjected to probit transformation (Zuk, 1989) prior to 
statistical analysis. 

All investigated traits, including gestation length, were successive calving (p w xl), calving season (p s xl), housing 
subjected to explorative analysis with the involvement of system (p w xl), calf sex (p P xl); 
multiple analysis of variance and the GLM procedure (SAS, 



The genetic parameters of gestation length and 
reproductive traits were estimated using Gibbs sampling 
method and the sire model. The components of 
(co-)variance were determined with the involvement of 
THRGIBBS1F90 (Tsuruta and Misztal, 2006) software and 
bivariate linear-threshold models. GL was an obligatory 
component in bivariate models, and the second trait was one 
of the investigated reproductive attributes: SB, CE and PE. 
The components of (co-)variance were estimated by Gibbs 
sampling to generate 50,000 samples, of which 10,000 were 
representative of the burn-in period. The number of rejected 
initial samples was computed with the use of 
POSTGIBBS1F90 (Tsuruta and Misztal, 2006) application. 
The results produced by neighbor samples were correlated 
(auto-correlated), and genetic parameters were determined 
based on data for every 100th sample. Standard errors of 
(co-)variance components were calculated as standard 
deviation of component values produced by successive 
samples. 

Other statistical sources of variability were determined 
in the GLM procedure, subject to the analyzed dependent 
variable, therefore, varied linear models were applied: 

GL: y = X w p w +X s p s +X P p P +Z s s+e 

SB: 1 = XuPu+XpPp+ZsS+e 

CE: 1 = XnpH+Xwpw+Xupu+Xppp+Z^+e 

PE: 1 = X H pH+X w pw+Xsps+X P p P +Z s s+e 

y = nxl observation vector; 1 = vector of liabilities; 
Ph, Pw, Ps, Pu, Pp = fixed effect vectors: herd (ph x 1), 
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s = qxl vector of random additive genetic sire effects; 

X H , X w , X s , Xu, X P = incidence matrices for fixed 
effects: herd (nxp H ), successive calving (nxp w ), calving 
season (nxp s ), housing system (nxp y ), calf sex (nxp P ); 

Z s = nxq incidence matrix for random additive genetic 
sire effects; 

e = nxl vector of random errors. 

It was assumed that the liabilities conditional on all of 
the effects were independent and normally distributed. 

The studied population's genetic parameters were 
determined using the following formula: 

Heritability h 2 = 4x °^ , 

genetic correlation cov s , 

V°s, XO s 2 

phenotypic correlation r cov s + cov e 

2 2 

where, o s = additive genetic sire variance, a E = 

error variance, where: cov s = genetic sire covariance. 

Genetic parameters were estimated in two replications. 
At the first stage, 47 half-sib groups were created based on 
the sires of calved cows. In the second step, 24 half-sib 
groups were created in view of the sires mated to cows. The 
minimum size of a half-sib group was 15 individuals. 

RESULTS AND DISCUSSION 

The average gestation length was 279.1 (±4.79) days 
(Table 1), and was shorter than determined for Holstein 
cattle by Nadarajah et al. (1989) 281.3 days, Silva et al. 
(1992) 280 days, and Norman et al. (2009) 281.6 days. In a 
study by Hansen et al. (2004), the average gestation length 



Table 1. Statistical characteristic for investigated traits 



Variable 


Mean 


SD 


Skewness 


Gestation lenght (d) 


279.54 


4.79 


-0.029 


Calving ease 


1.60 


0.88 


1.36 


Stillbirths 


1.06 


0.24 


3.61 


Placental expulsion 


1.47 


0.72 


1.17 



in Danish Holstein cows was 278.5 days. Johanson and 
Berger (2003) reported average GL values in American HF 
cattle at 277.9 days. The gestation length of Guernsey, 
Holstein and Jersey cows in several Florida-based farms 
increased by four days on average over a period of 50 years. 
According to the authors Silva et al. (1992), this increase 
could be attributed to a steep rise in production levels. 
Shorter gestation was noted in the summer months, and 
these findings are consistent with the results of previous 
research which validated the effect of higher temperature 
and dietary changes on shorter GL (Hansen et al., 2004; 
Przysucha and Grodzki, 2009). In this study, a longer 
gestation period was observed in respect of older cows 
(Table 3). There was a significant difference between the 
first and the third calving, at one day on average. In a study 
analyzing only primiparous cows, Hansen et al. (2004) 
reported a difference of 3.5 days between heifers calved at 
the age of 22 and 36 months. In this study, the increase in 
calf weight significantly contributed to a rise in GL values. 
In the group of the heaviest calves (>40 kg), GL was 3.9 
days longer on average in comparison with the group of the 
lightest calves (<34 kg). Nadarajah et al. (1989) reported a 
difference of 3.7 days between cows with large and small 
fetuses. In this study, GL values were also significantly 
affected by the sex of the fetus, and they were 1.8 days 
longer for male fetuses. According to Silva et al. (1992) and 
Hansen et al. (2004), GL was 1.1 days longer for male 
calves. 

Gestation length increased with the degree of calving 
difficulty (Table 4). Difficult and very difficult calvings 
were characterized by a significant increase in GL values. A 



Table 2. Distribution (%) of calving ease, stillbirths and placental expulsion in respect of investigated factors 



Calving ease Stillbirths Placental expulsion 







1 


2 


3 


4 


1 


2 


1 


2 


3 


Calving season 


Fall-winter 


66.67 


21.11 


8.89 


3.33 


95.33 


4.67 


70.44 


19.11 


10.44 




Winter-spirng 


55.51 


24.25 


13.43 


6.81 


92.99 


7.01 


62.53 


22.04 


15.43 




Summer 


62.66 


24.03 


7.14 


6.17 


92.53 


7.47 


64.29 


22.73 


12.99 


Successive calving 


First 


45.25 


30.28 


15.85 


8.63 


91.90 


8.10 


66.73 


25.00 


8.27 




Second 


68.88 


20.82 


7.09 


3.20 


94.51 


5.49 


65.90 


17.39 


16.70 




Third 


84.13 


10.71 


3.17 


1.98 


96.43 


3.57 


63.49 


19.05 


17.46 


Calf sex 


1 (heifer calves) 


63.34 


23.99 


8.95 


3.72 


95.61 


4.39 


70.10 


21.28 


8.61 




2 (bull calves) 


59.40 


22.26 


11.43 


6.92 


92.03 


7.97 


61.95 


21.05 


16.99 


Husbandry system 


1 (tie-stall) 


54.82 


26.76 


11.37 


7.05 


92.23 


7.77 


70.36 


19.57 


10.07 




2 (free-stall) 


69.22 


18.51 


8.90 


3.38 


95.55 


4.45 


60.14 


23.13 


16.73 
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Table 3. Least square means (LSM) of gestation length subject to 
calving season, calving number, calf weight and sex and housing 
system 

Gestation length (d) 





n 






Calving season 








Fall-WIIILCI 


yvjvj 


970 1 




(1 October-31 January) 








Winter-spring 


998 


279.3 


4.84 


(1 February-31 May) 








Summer 


616 


278.5 


4.88 


(1 June-30 September) 








Successive calving 








First 


1,136 


278.6 a 


4.89 


Second 


874 


279.3 


4.74 


Third 


504 


279.6 b 


4.55 


Calf weight (kg) 








<34 


626 


278.0 A 


4.37 


35-37 


670 


279.2 B 


4.28 


38-40 


558 


280.7 C 


4.60 


>40 


300 


281.9 D 


4.72 


Calf sex 








Heifer calves 


1,184 


278.1 A 


4.39 


Bull calves 


1,330 


279.9 B 


4.97 


Housing system 








Tie-stall 


1,390 


279.0 


4.72 


Free-stall 


1,124 


279.1 


4.87 



Values followed by different 
significantly: capital letters p<0. 



letters in columns within factors differ 
01, small letters p<0.05. 



significant correlation between GL and CE was 
demonstrated. Difficult and very difficult calvings that 
required veterinary assistance had a 15.67% share of all 
births. The highest share of difficult and very difficult 
calvings was reported in the group where GL exceeded 283 
days, whereas the easiest calvings were observed in the 
group of cows whose GL ranged between 275 and 277 days. 
In a study by Cervantes et al. (2009), the GL of cows 
experiencing calving difficulties was 2.1 longer in 
comparison with cows that calved without assistance. 
Philipsson (1976) demonstrated a correlation between 
difficult calving and higher GL values in Swedish Friesian 



heifers. 

Both prolonged and shortened gestation contributed to a 
significant increase in the number of stillbirths (Figure 1). 
The lowest SB rate (3.23%) was observed in respect of GL 
values in the range of 275-277 days, whereas the share of 
calves that were born dead and died within 24 h after birth 
increased three-fold (to 10.48%) when GL exceeded 283 
days. Johanson and Berger (2003) and Philipsson (1976) 
reported a non-linear correlation between GL and SB. In 
this study, optimal gestation length in Holstein cows was 
determined in the range of 275-277 days based on calving 
ease and stillbirth values. In a study evaluating more than 
4,000 Friesian cows characterized by average gestation 
length of 278.7 days, Johanson and Berger (2003) found 
that the optimal GL that minimized the risk of SB was 282 
days. Deviations from the optimal GL value led to a 
significant increase in SB rates. According to Norman et al. 
(2011), average GL values optimized lifetime productivity, 
calving ease, stillbirth rates and calving-to-first-service 
interval. 

In the discussed experiment, genetic parameters were 
estimated with the use of Gibbs sampling and bivariate 
linear-threshold models. A similar approach was adopted by 
Lee et al. (2002), Ramirez-Valverde et al. (2001), Varona et 
al. (1999) who investigated the adequacy of threshold-linear 
models for estimating breeding value based on calving ease, 
birth weight and gestation length. The above authors 
concluded that linear-linear models and threshold-linear 
models produced a similar sequence of breeding values for 
birth weight and gestation length. A lower correlation was 
determined with regard to calving ease (threshold trait). 
Threshold-linear models delivered more accurate 
estimations than linear- linear models. 

The heritability of gestation length that accounted for 
direct sire of fetus effects varied in the range of 0.201-0.210 
(Table 6). Heritability values determined in view of indirect 
maternal effects were visibly lower at 0.054-0.073 (Table 5). 
GL heritability indicators reported in this study were 
generally lower than the values reported by other authors 
(Hansen et al., 2004; Crews, 2006; Norman et al., 2009; 
Johanson et al., 2011). 

The results of this study as well as the cited authors' 
findings were characterized by significant variance which 



Table 4. Calving ease vs. gestation length (x 2 = 35.1; p<0.01) 



Gestation length (d) 



Calving ease 




<275 


275-277 


278-280 




281-283 




>283 


n 


% 


n 


% 


n 


% 


n 


% 


n 


% 


Unassisted 


254 A 


55.95 


408 B 


73.12 


344 B 


66.70 


328 B 62.12 


234 A 


51.09 


Little intervention 


134 a 


29.52 


88 b 


15.77 


90 b 


22.87 


122' 


23.11 


118 a 


25.76 


Difficult 


44 


9.69 


38 


6.81 


54 


10.47 


46 


8.71 


76 


16.59 


Very difficult 


22 


4.85 


24 


4.30 


28 


5.43 


32 


6.06 


30 


6.55 



Values followed by different letters in rows differ significantly: capital letters p<0.01; small letters p<0.05. 
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Figure 1. Gestation length vs. perinatal mortality and unassisted calving. 



can be attributed to numerous factors, including the 
classification model (linear or threshold), type of random 
effects, influence of breed or successive lactation (Hansen 
et al., 2004; Crews 2006; Cervantes et al., 2009; Norman et 
al., 2009; Johanson et al., 2011). 

The type of random effects in the classification model 
had a significant influence on the estimated heritability 
indicators as regards calving ease (Tables 5 and 6). When 
the random effect of sires of calved cows was taken into 
account, stillbirth heritability increased two-fold, whereas 
the heritability of placental expulsion was reduced more 
than three-fold, and the heritability of calving ease 
decreased more than two-fold. 

The model accounting for the random effect of sires of 
calved cows demonstrates moderate positive genetic 
correlations between GL vs. SB and CE (Table 5). The 
genetic optimization of gestation length could have a 



negative influence on the breeding value of bulls by 
increasing stillbirth frequency and calving difficulty. The 
use of a classification model accounting for the sire of calf 
model produced correlations between GL vs. stillbirth rates 
and calving ease that were characterized by an identical 
pattern but lower strength (Table 6). Similar genetic 
correlations (direct effect) between gestation length vs. 
stillbirth rates (0.18) and calving difficulty (0.38) were 
reported by Hansen et al. (2004) with the use of a linear- 
threshold model. Subject to the type of random effect, 
genetic correlations with a different sign were reported for 
correlations between GL and placental expulsion. A 
negative correlation was determined for the random effect 
of sires of calved cows, while a positive correlation was 
reported for the effect of sires of live-born calves. 

The coefficients of all phenotypic correlations between 
gestation length and calving ease were positive but low. The 



Table 5. Heritability, genetic and phenotypic (SE) correlation between gestation length and calving ease subject to indirect sire effects 



Parameter Gestation length 


Stillbirths 


Placental expulsion 


Calving ease 


Heritability 0.054-0.073 (0.032) 
Genetic correlation 
Phenotypic correlation 


0.420 (0.157) 
0.555 (0.321) 
0.088 (0.051) 


0.155 (0.067) 
-0.340 (0.415) 
0.040 (0.025) 


0.211 (0.085) 
0.625 (0.294) 
0.108 (0.030) 


() standard deviation is shown in brackets. 








Table 6. Heritability, genetic and phenotypic (SE) correlation between gestation length and calving ease subject to direct sire effects 


Parameter Gestation length 


Stillbirths 


Placental expulsion 


Calving ease 


Heritability 0.201-0.210 (0.062) 
Genetic correlation 
Phenotypic correlation 


0.195 (0.099) 
0.179 (0.453) 
0.078 (0.044) 


0.560 (0.182) 
0.149 (0.244) 
0.046 (0.029) 


0.552 (0.169) 
0.242 (0.161) 
0.113 (0.032) 



() standard deviation is shown in brackets. 
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above suggests that prolonged gestation resulted in greater 
calving difficulty, higher stillbirth rates and longer placenta 
expulsion time. 

An analysis of standard deviations of the investigated 
parameters revealed relatively high values (Tables 5 and 6). 
For this reason, modifications of gestation length should be 
introduced with great caution during attempts to improve 
calving ease and reduce stillbirth rates. 

CONCLUSIONS 

The following factors contributed to prolonged 
gestation in Holstein-Friesian cows: increasing cow's age, 
male fetuses and growing fetus weight. Optimal gestation 
length was determined in the range of 275-277 days based 
on calving ease and stillbirth values. In view of direct sire 
of fetus effects, GL heritability was found to be four-fold 
higher in comparison with heritability values produced by 
indirect maternal effects, the model accounting for the 
random effect of sires of calved cows demonstrated 
moderate positive genetic correlations between GL vs. SB 
and CE. The above suggests that genetic optimization of 
gestation length could have a negative influence on the 
breeding value of bulls by increasing stillbirth frequency 
and calving difficulty. The use of a classification model 
accounting for the sire of calf model produced correlations 
between GL vs. SB and CE that were characterized by an 
identical pattern but lower strength. The investigated 
parameters were marked by relatively high standard 
deviations, implying that GL modifications should be 
introduced with great caution during attempts to improve 
CE and SB values. 
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